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ABSTRACT 
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density  profile  was  Gaussian  shaped.  The  second  method  used  the  Abel  Integral  Trans¬ 
formation,  requiring  no  pre-conditions  on  the  density  profile  except  that  it  was 
eylindrically  symmetric.  While  used  extensively  in  plasma  spectroscopy,  this  technique 
is  not  generally  well  known  for  analyzing  interferometric  data.  The  results  show  that  for 
a  fixed  laser  energy  long,  shallow  reduced-density  channels  were  formed  in  gas  mixtures 
of  low  ammonia  concentration  and  short,  deep  channels  were  formed  in  gas  mixtures 
of  high  ammonia  concentration.  These  results  qualitatively  agree  with  an  earlier  exper¬ 
iment  in  which  gas  mixtures  of  nitrogen  and  sulfur-hexaflouride  were  used.  Both  exper¬ 
iments  support  the  concept  of  reduced-density  channel  formation  and  lay  the  foundation 
for  future  studies  of  relativistic  electron  beam  propagation  for  application  in  the  Strate¬ 
gic  Defense  Initiative. 
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I.  IMRODL  C  I  ION 


The  Strategic  Defense  Initiative  •  sDl  •  „ 


dear  weapons  obsolete.  One  of  the  mans  areas  be  mg 
goal  is  the  use  of  intense  high  cncrgv  electron  beams 
these  beams,  it  was  found  that  interaction  w.th  the  at:: 
propagation.  To  improve  beam  propagation,  a  met  ho 
mospheric  density  along  the  path  of  the  beam  1  h:s  met 
formation,  consisted  of  using  a  laser  to  pre-vond.t.  tt 
gion  of  reduced  demits  along  the  intended  path  c:  t:.c  c 
Studies  of  Reduced  Demit*.  Junnfs  ••'•etc  .  . 


the  Strategic  Defense  Initiative  [Rets  ■  1 


Alamos  National  Laboratory.  to  conduct  experiments  m 
of  experiments  were  conducted  m  '.ditch  variable  lact  c: 
incident  on  target  cells  containing  gas  mixtures  One  g: 
cells  containing  mixtures  of  nitroecn  and  si.ifar-hcx.if  ■ 
second  group  of  experiments  which  used  target  <<  b  m 
trosren  and  ammonia  are  di'cU"c\;  here  \  'J.cn  -t:> 


Figure  1 


A  <'02  laser  beam  incident  >ui  a  tarcet  >e  '. 


creates  a  re  cion  of  reduced  demit*,  when  the  Mi’-  rut;.. 


characteristics 


the  rediiccd-dcnsitv  ciiunne.  a: 


interferoerams  taken  alter  the  Li'-or  was  fired 


The  (02  laser  used  m  exnerm.ent  w as  fur.a'*'c 


3IM>  joules.  The  ho  r  pul  sc  s  exhibited  t  time  p" 
guassian  sp  ina;  pr.  :.,c  \ppm\;m  itch.  ;;  pc.c 


tlie  aperture  and  enorg-ueJ  the  gas  m  the  ceh.  that  -.s 
available  lor  channel  formation.  As  the  '"cam  tr.nc:« 


mixture  ol  Nile  and  N2.  a  recion  cf  reduced  der.s;* 


i,rce  intcrteromcters  were 


t  See  f  i  .ure  I  >. 
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Pictures  o t  the  .uterteren.c 


the  laser  beam  cas  t  t  .  ■ .  '  __  .  .  .  „  .  -  . 

t.)  detern  :m :  'he  d:  oth  arid  .  ..  me  _  .  .  ....  I  ■  u.  .  .  ...  • 

to  determine  me  .hara.tcris'-es  o!  dcm.i .  ..... 

I  tie  first  technique,  trie  (lacsiu::  ::.c‘..'u r  ass..:r.cd  me  dcr. ».'.■•  pro.de 
from  laser  absorption  had  a  t»au>s;art  s.’.a  <<  L.c  sc.  cud  tc  rm...  ..i  „cu..-;\-c  : .  - 
rtune  the  channel  characteristics  rcqu.rcd  ::.c  a>*  ..rr.pt.v-r.  :.  the  he  mt\  pr . 
cvhndricalh  s-.mmctrio  the  \b;.  .'.terra,  tr-m.  rural..,:.  _  ac..  cm  s ;m.m... 
no stic  technique.  ‘.vax  adapted  tor  use  •.*. rr.tertcr .  me'.?.*  d..ra  1  h.s  ..«.’.»pt..i:.  . 

hesed  to  he  one  at  the  f.rst  apph-at.  urn  the  \  •<:  m  .j.tcrUti  n 


:i.s  thesis  details  the  pr  m;  :  . 


'eristics  from  intcricrcnce  p  i"  >  rc.  rded  the  c\pc.*.::.e:.t  arid  .  mm..' 

thorough  analysis  of  all  the  d-.h.  an  the  *c-.  om  :•  up  !  ■  :.c  !  t  s  \..;u..s  espcr.m 
I  hose  results,  described  m  Chapter  v.  .;;  md.catc 


[  or  a  fixed  mere'-  lc. c!,  increasing  h-c  >. on  c: .t r oj  \H  :  result*  » 
shorter  reduecd-den  .its  channels 


1  ucepc: 


•  [  or  a  (i.iu"i  m  '!: aped  m.terter.me  p  •  .. 

'•ides  a  quick  and  .u, urate  n  ear.'  :  ..."cm..:.  :  c  *. i.c  reduce  *• 
characteristics 


. . c  r c ».* •  * v. c  ■hUi’1.'’  k .* - •  j 


These  results  ire  supported  he  the  N2-M  *•  experimental  result'  and  a  hi  he  mid  to  s 
the  propagation  of  intense  electron  beams  throuch  ti.c  atm.'  sphere. 


II.  DESCRIPTION  OF  THE  REDIC  ED- DENSITY  CHANNEL 

EXPER1MEN  I 


A.  EQUIPMENT 


A  schematic  of  the  experiment  is  shown  in  figure  1.  A  tunable  electron  beam  con¬ 
trolled  CQ2  laser  with  an  llxll\l')0  cm’  active  volume  produced  pulses  of  I'M'  uni  ra¬ 
diation  at  energies  of  up  to  300  joules  hasing  a  3.3  us  pulse  width.  Cpon  leasing  the 
laser,  the  radiation  experienced  multiple  reflections  as  it  passed  through  a  nitrogen  filled 
cylindrical  tube  10  m  long.  The  reflections  served  to  increase  the  optical  path,  hence  the 
laser  beam  traveled  approximately  33  m  before  exiting  the  tube  and  entering  the  test  cell. 
1  he  increased  length  aided  in  focusing  the  beam:  the  nitrogen  in  the  cylinder  presented 


the  beam  from  being  scattered  by  particulate  matter.  Since  nitrogen  does  not  absorb  I't.o 
win  radiation  the  laser  beam  approached  the  apperture  without  any  significant  degrada¬ 
tion.  One-third  of  the  initial  laser  energy  passed  through  the  aperture  and  was  incident 
upon  the  gas  at  the  entrance  of  the  target  cell.  A  ty  pical  60  joule  pulse  entering  the  target 
cell  had  an  estimated  peak  on-axis  llucncc  of  5*>  Joules  per  an2.  Entering  the  cell  with 
a  full  width  at  half  maximum  of  ~10  mm.  the  beam  was  aligned  so  that  the  peak  of  its 
Gaussian  profile  was  along  the  longitudinal  axis  of  the  cylindrical  tost  cell.  (Ref.  S:  pp. 
4-3], 


The  test  cell,  a  cylinder  36  cm  in  diameter  and  230  cm  long,  contained  a  mixture  of 
N1I3  and  N2.  Experiments  in  the  test  cell,  were  conducted  on  the  following  mixtures: 

•  5  torr  NII3  and  305  torr  N2 


•  15  torr  NII3  and  5S5  torr  N2 

•  50  torr  NII3  and  550  torr  N2 


•  200  torr  N1I3  and  dOOtorr  N2 


The  test  cell  was  evacuated  before  each  experiment.  Subsequently  it  was  filled  to  the 
desired  pressure  of  3N 1 13  and  then  brought  to  atmospheric  pressure  (600  torr  for  Eos 
Alamos*  with  pressurized  N2. 

The  XI 13  molecules  absorb  10.6  /an  radiation  (Refs.  9.  1(>  and  1 1].  As  the  beam  il¬ 
luminated  the  cell,  the  absorbed  energy  was  transferred  to  the  other  molecules  of  the  gas 
mixture  tlv-migh  collisions.  1  he  adiabatic  expansion  of  these  excited  molecules  produced 
a  reduced  density  channel  which  was  centered  on  the  longitudinal  axis  of  the  test  ceil. 


3 


An  in  depth  description  of  the  adiabatic  expansion  process  is  developed  in  reference 
S  pages  5  through  10.  Included  in  that  reference  are  calculations  which  determine  the 
required  absorbed  lluence  to  produce  a  unform  half  density  channel  2  m  in  length.  While 
the  values  obtained  therein  apply  to  600  torr  N2.  they  apply  to  the  mixtures  used  in  this 
experiment  due  to  the  prominence  of  the  N2  contained  in  each  gas  mixture.  In  short  the 
tluences  provided  by  the  C02  laser  were  sufficient  to  produce  2  m  long  channels  which 
were  1  cm  wide  and  had  0.1  times  the  ambient  density. 

In  order  to  observe  this  reduction  in  density,  windows  were  installed  along  the  walls 
of  the  cell.  These  windows  permitted  the  use  of  optical  interferometers,  one  of  whose 
beams  traversed  the  test  cell.  The  density  profile  of  the  gas  mixture  at  different  locations 
was  then  determined  by  analyzing  the  interference  patterns. 


B.  INTERFEROMETERS 

Three  interferometers  were  installed  at  different  positions  along  the  test  cell  (See 
Figure  1).  These  Mach-Zehnder  type  of  interferometers  split  the  Ile-Ne  laser  (633nm) 
into  two  separate  beams.  One  beam  traversed  across  the  path  of  the  C02  laser  beam 
and  the  other  beam  traversed  around  it.  Figure  3  is  a  photograph  of  the  pattern  resulting 
from  interference  of  both  beams  prior  to  the  energizing  of  the  C02  laser  and  subsequent 
channel  development.  This  pattern  is  due  to  non-alignment  of  the  mirrors  used  in  the 
Interferometer  [Ref.  12).  Filling  the  test  cell  with  the  \II3  and  \2  mixture  and  energizing 
the  C02  laser  developed  a  phase  shift  between  the  two  beams.  This  phase  shift  was  due 
to  a  decrease  in  the  optical  path  length  of  the  beam  passing  across  the  path  of  the  (.'02 
laser.  Figure  4  is  a  photograph  of  the  resulting  pattern. 

Photographs  were  taken  of  each  interference  pattern  observed.  These  photographs, 
interferocrams.  were  processed  in  order  to  determine  the  density  profile  The  processing 
required  photographic  enlargment.  digitizing  and  numerical  anahsiv  fhe  following  sec¬ 
tion  is  devoted  to  describing  terms  used  in  anaix/ing  the  interferocrams. 

C.  INTERFEROCRAMS 


An  mtcrlerogram  is  a  photographic  record  of  a  pattern  produced  1". 
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Top  and  Bottom  of  interferogram. 

The  edge  of  the  intcrferograrn  towards  which  the  peaks  of  the  fringe  line  are 
displaced  is  the  top.  The  bottom  of  the  interferogram  is  the  edge  directly  opposite 
the  top. 

2  cm  reference  mark. 

In  order  to  establish  a  scale  between  the  interferograms.  a  2  centimeter  refer¬ 
ence  mark  was  devised  and  located  in  the  field  of  view  of  the  camera.  This  device, 
a  periodic  2  cm  square-tooth  pattern  was  used  for  each  interferogram.  Portions  of 
the  device  are  visible  in  Figures  3  and  4.  In  Figure  5  the  device  is  shown  at  the  top 
of  the  sketch.  The  distance  from  the  right-most  corner  of  one  tooth,  to  the  right¬ 
most  comer  of  the  next  tooth  was  2  cm. 

Fringe  line: 

The  interference  pattern  showns  in  Figures  3  and  4.  also  sketched  in  Figure  5, 
consist  of  alternating  regions  of  black  and  white.  Although  either  region  could  be 
used  to  define  a  fringe  line,  the  black  region  is  more  clearly  defined  visually.  Thus 
throughout  this  description,  a  fringe  line  refers  to  a  line  in  the  black  region  which 
is  halfway  between  two  white  regions.  Figure  5  shows  a  sketch  of  two  typical  fringe 
lines. 

Fringe  number 

The  fringe  number  was  arbitrarily  chosen  to  be  the  number  of  fringe  lines  from 
the  bottom  of  the  interferogram.  In  Figure  4,  the  third  fringe  is  found  by  starting 
at  the  first  black  region  at  the  bottom  of  the  interferogram.  The  fringe  line  in  this 
region  is  fringe  number  one.  The  next  black  region  encountered  as  one  moves  from 
bottom  to  top  is  fringe  number  two.  contains  fringe  number  two.  The  next  region 
contains  fringe  number  three  and  so  on. 

Zero  fringe  line. 

A  local  term  used  to  denote  a  straight  line  passing  from  the  left-most  edge  of 
a  fringe  line  to  the  right-most  end  of  the  same  fringe  line.  This  line  corresponds  to 
the  interference  pattern  which  would  be  present  if  there  were  no  fringe  displacement 
(i.e.,  no  absorbing  gas  in  the  test  cell).  Two  of  these  zero  fringe  lines  are  annotated 
in  Figure  5. 

Fringe  displacement. 

The  distance  measured  in  fringe  lines  above  the  zero  fringe  line.  It  is  due  to 
the  decreased  optical  path  length  created  when  the  gas  mixture  absorbs  10.6  fim 
radiation.  To  determine  fringe  displacement,  one  must  find  the  distance  between 
successive  zero  fringe  lines.  A  ratio  of  these  numbers  yields  a  scale  factor  from 
which  any  distance  measured  along  the  vertical  can  be  converted  to  a  fringe  value. 
In  Figure  5.  fringe  number  one  has  zero  fringe  displacement  at  the  far  left  and  the 
far  right.  The  peak  of  this  fringe  line  is  approximately  1.25  fringe  lines. 

Vertical  reference  line. 

A  vertical  reference  line  was  drawn  connecting  the  peaks  of  the  fringe  lines. 
This  line  was  parallel  to  the  vertical  faces  of  the  square  tooth  but  was  not  neces¬ 
sarily  perpendicular  to  the  zero  fringe  line. 


D.  summary  of  experiment 

In  this  chapter  the  experiment  and  its  recorded  results  are  described  I  he  ne\ 
in  trying  to  determine  the  density  profile  and  subsequent!;.  the  chaiacteii't^v 
reduced-density  channel,  is  to  interpret  the  imerl'erograms  Tins  interpret.!:. on  : 
ducted  in  the  next  chapter. 


III.  DATA  REDUCTION 


Photographs  of  the  interference  patterns  were  taken  shortly  after  the  CO 2  Laser 
energized  each  gas  mixture.  These  interferograms  presided  a  record  from  which  the 
fringe  data,  the  relationship  between  fringe  displacement  and  displacement  along  the 
zero  fringe  line,  would  be  obtained  and  the  on-axis  density  determined.  Once  the  fringe 
data  was  obtained  the  density  profile  and  subsequently  the  reduced-density  channel 
characteristics  could  be  determined.  This  chapter  describes  how  the  fringe  data  was  ex¬ 
tracted  from  the  interferogram  and  placed  into  data  files.  These  files  will  he  analyzed 
later  by  the  Gaussian  and  Abel  techniques  to  determine  the  characteristics  of  the 
reduced-density  channel  formed  in  each  gas  mixture. 

A.  PHOTOGRAPHIC  ENLARGEMENT 

Determining  the  fringe  data  required  making  measurements  directly  on  the 
interferogram.  To  facilitate  this  procedure,  each  original  interferogram.  a  photograph 
measuring  approximately  3  cm  by  3  cm.  was  enlarged.  The  enlargement  process  was 
accomplished  using  a  specially  modified  35  mm  camera  fitted  with  a  zoom  lens.  The 
image  observed  in  the  field  of  view  of  the  camera  was  converted  to  an  electronic  signal 
and  simultaneously  passed  to  a  video  screen  and  Tektronix  photo-copier.  While  moni¬ 
toring  the  video  screen,  the  focus  of  the  camera  was  adjusted  until  a  clear  image  of  the 
2  cm  reference  marks  and  interference  pattern  filled  the  video  screen.  The  copy'  button 
of  the  photo-copier  was  then  depressed,  resulting  in  an  enlarged  photograph  of  the 
interferogram.  Some  interferograms  were  not  capable  of  being  enlarged  with  the  2  cm 
reference  marks  visible.  For  those  interferograms.  a  scribe  was  used  to  mark  the  refer¬ 
ence  points  on  the  orignal  photograph.  The  scribe  marks  enhanced  the  image  of  the 
reference  marks  and  made  them  visible  in  the  enlarged  interferogram. 

Each  enlarged  photograph  therefore  contained  a  set  of  marks  which  measured  2  cm 
in  the  original  interferogram.  The  distance  between  these  marks  varied  with  each 
interferogram  as  a  result  of  the  focusing  adjustments  made  during  the  enlargement 
process.  Thus,  all  distances  could  be  scaled  to  the  common  2  cm  reference  mark  and 
continuity  was  maintained  between  photographs. 


Having  enlarged  each  interferogram.  the  next  step  in  determining  the  density  profile 
and  subsequently  the  characteristics  of  the  reduced-density  channel  was  to  digitize  the 
enlarged  photographs.  Before  digitizing  the  interferograms.  a  one-time  check  of  the 
combined  effects  of  digitizing  and  enlarging  was  conducted  by  comparing  an  original 
drawing  to  an  enlarged  copy  of  the  same  drawing.  A  line  was  drawn  on  a  piece  of  graph 
paper  and  eight  easily  identified  points  were  then  digitized  (See  Figure  6).  The  drawing 
was  then  enlarged,  using  the  process  outlined  above,  and  the  same  eight  points  on  the 
enlargement  were  digitized.  The  original  and  resulting  coordinates  are  listed  in  Table  1. 
A  review  of  this  Table  illustrates  that  the  maximum  difference  between  the  original  and 
the  enlarged  copy  are  within  .04  cm  of  each  other.  The  mean  difference  between  the  or¬ 
iginal  and  the  enlarged  copy  was  no  greater  than  .02  cm  for  both  the  vertical  and  hori¬ 
zontal  components.  It  will  be  shown  in  chapter  6  that  the  errors  generated  by  the 
enlargement  process  did  not  significantly  effect  the  results.  The  next  step  in  determining 
the  density  profile  was  to  prepare  the  interferograms  for  digitization. 

B.  DIGITIZING  PREPARATIONS 

In  preparation  for  digitization,  each  enlarged  interferogram  was  reviewed  and  the 
following  items  identified  by  marking  on  the  photograph.  Figure  7  is  a  sketch  of  an 
inteferogram  prepared  for  digitizing. 

•  2  cm.  reference  marks. 

The  2  cm  reference  mark  was  identified  by  drawing  a  line  from  the  lower  right 
corner  of  the  left-most  square  tooth.  Rl.  to  the  same  comer  of  the  next  square 
tooth.  R2. 

•  Fringe  lines  to  be  digitized. 

Two  fringe  lines  per  interferogram  were  chosen  to  be  digitized.  These  fringe 
lines  were  normally  picked  based  on  their  clarity  and  proximity  to  the  center  of  the 
interferogram.  Three  fringe  lines  labelled  FI.  F2.  and  F3  are  shown  in  Figure  7. 

•  Vertical  reference  line. 

A  vertical  reference  line  was  drawn  connecting  the  peaks  of  the  fringe  lines. 
This  line  was  parallel  to  the  vertical  faces  of  the  square  tooth  but  was  not  neces¬ 
sarily  perpendicular  to  the  zero  fringe  line. 

•  Zero  fringe  lines. 

Zero  fringe  lines  were  drawn  for  many  of  the  fringe  lines  of  each  interferogram. 
These  lines  provided  the  basis  for  determining  the  fringe  displacement  and  subse¬ 
quently  the  on-axis  density.  Three  zero  fringe  lines  have  been  drawn  in  Figure  7. 

Having  prepared  each  interferogram  for  digitizing,  the  next  step  was  translating  the 
photographic  information  available  from  the  interferogram  to  a  form  which  would  allow 
numerical  processing. 


C.  DIGITIZATION 


Converting  the  photographic  image  of  the  interference  pattern  into  a  form  which 
allows  mathematical  computation  is  the  process  of  digitizing. 

The  first  step  in  the  process  was  to  establish  a  rectangular  coordinate  system  for  the 
interferogram.  This  was  accomplished  by  designating  the  left-most  point  on  the  fringe 
line  as  the  origin.  This  point  is  labelled  P(l)  in  Figure  7.  The  vertical  axis  was  made 
parallel  to  the  vertical  reference  line  of  the  interferogram  and  the  horizontal  axis  was 
chosen  perpendicular  to  this.  With  this  reference  system  established  any  point  on  the 
interferogram  could  be  designated  by  a  vertical  and  horizontal  component.  The 
HP-91 11A  Graphics  Tablet  was  used  to  assign  numerical  values  to  these  components. 
The  numerical  values  are  in  digitizing  units,  where  1  digitizing  unit  is  equal  to  0.025 
millimeter  [Ref.  13].  The  Graphics  Tablet  can  be  thought  of  as  an  electonic  grid  with  its 
own  rectangular  reference  system  and  units.  The  wand  used  with  the  Tablet  had  a  spring 
loaded  tip.  With  the  tip  in  contact  with  the  Tablet  and  the  spring  compressed,  the  posi¬ 
tion  of  the  wand  relative  to  the  reference  system  of  the  Graphics  Tablet  was  registered. 
This  information  was  then  passed  to  the  HP-S5  computer  and  stored  in  a  data  file. 

Securely  fastened  to  the  Tablet  using  adhesive  tape,  the  interferogram  was  ready  to 
be  processed.  This  consisted  of  digitizing  the  following  characteristics  of  the 
interferogram. 

•  2  cm  reference  mark. 

Digitizing  the  right-most  corner  of  one  square  tooth  and  then  the  same  corner 
of  the  next  tooth  allowed  the  distance,  in  digitizing  units,  to  be  determined.  This 
distance  corresponded  to  2  cm  on  the  interferogram.  Refer  to  points  labelled  R1 
and  R2  on  Figure  7. 

•  Distance  between  zero  fringe  lines. 

Digitizing  the  intersection  of  the  vertical  reference  line  and  a  zero  fringe  line 
and  the  vertical  reference  line  with  different  zero  fringe  line  allowed  a  scale  factor 
to  be  determined.  The  distance  between  these  two  digitized  points  represented  the 
number  of  digitizing  units  per  fringe.  Refer  to  points  R3.  R4  and  R5  of  Figure  7. 
Having  previously  determined  the  number  of  digitizing  units  per  2  cm.  a  value  for 
the  number  of  fringes  per  cm  was  determined. 

•  Fringe  line. 

All  fringe  lines,  two  per  interferogram.  were  digitized  in  the  same  manner. 
Starting  at  the  left-most  point  of  the  fringe  line,  point  P(  1 1  Figure  7.  and  proceeding 
towards  the  right,  approximately  35  points  per  line  were  digitized.  (See  fringe  line 
FI  of  figure  7.)  flic  last  point  digitized  was  the  right-most  point  of  the  fringe  line. 
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along  the  zero  fringe  line,  some  data  manipulation  was  needed  Mas  manipulation 
consisted  of  correcting  the  digitized  tungc  ime  tor  the  angular  Jiil'cieme  between  the 
zero  fringe  line  and  the  horizontal  reference  line.  At  tins  punt  it  is  important  to  recall 
that  the  fringe  \aluc  is  determined  h\  the  distance  the  fringe  line  is  awa>  from  the  zero 
fringe  linc-not  the  horizontal  reference  hnc  1  iic  toil  awing  o-ut!..,es  ti.e  method  used  m 
writing  a  computer  program  which  con'. cried  digitized  data  ;nt  >  lunge  data. 

D  FRINGE  data 

Ihe  fringe  displacement  as  a  function  of  distance  along  the  zero  :> m.e  hue  defines 
a  rectangular  coordinate  sv'tcm.  In  order  to  obtain  tins  relationship,  a  cmndai  ite 
transformation  was  made,  lhts  transformation  required  aligning  tiie  itouzontai  a.\;s  of 
the  digitized  coordinate  ssstem  with  the  zero  fringe  line.  The  process  outlined  below 
follows  closely  the  development  of  a  computer  program  which  performed  the  transfor¬ 
mation.  To  assist  in  this  development,  the  following  terms  arc  defined. 

•  J  is  a  number  assigned  to  the  point  digitized.  J  =  1  lor  the  first  point  digitized.  .1  --  2 
for  the  second,  and  so  on. 

•  F  is  the  number  of  zero  fringe  lines  between  the  fringe  lines  digitized  for  hinge 
distance  calibration.  In  Figure  “.  the  number  is  1  if  R3  and  R4  were  digitized.  1 
equals  2  if  R3  and  Ro  were  digitized. 

•  V6  is  the  angle,  in  radians,  between  the  zero  fringe  line  and  the  horizontal  axis. 

•  D 1  is  the  distance,  in  digitizing  units,  between  the  2  cm  reference  marks. 

•  D2  is  the  distance,  in  digitizing  units,  between  zero  fringe  lines  digitized  for  a  ref¬ 
erence. 

•  Nl  is  the  last  point  of  the  fringe  line  digitized;  the  right-most  point. 

•  XI  is  the  horizontal  component  of  a  point  in  digitizing  units. 

•  Y1  is  the  vertical  component  of  a  point  in  digitizing  units. 

•  AtJ.l )  is  the  horizontal  component  of  a  digitized  point  J. 

•  A(J.2)  is  the  vertical  component  of  a  digitized  point  J. 

In  Figure  S.  one  of  the  many  fringe  lines  of  an  interferogram  has  been  sketched  and 
annotated  to  show  the  ecometrv  of  the  coordinate  transformation,  llavine  identified  the 
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geometric  relationships  in  the  interferogram.  the  next  step  in  generating  fringe  data  was 
referencing  all  digitized  points  to  the  origin.  This  referencing  was  accomplished  by 
subtracting  the  coordinates  of  the  origin  from  each  point  digitized.  As  shown: 

AT  =  .1(7.1)  -  .1(1.1)  i'.Ii 

Y\  =  A{J,2)  -  .1(1.2)  (3.2) 

The  next  step  was  to  determine  the  angle  between  the  zero  fringe  line  and  the  hori¬ 
zontal  axis.  This  was  accomplished  by  using  the  first  and  last  points  of  the  digitized 
fringe  line.  The  ratio  of  the  difference  of  vertical  components  to  horizontal  components 
was  the  tangent  of  the  desired  angle,  V6. 

A'6  =  /l(.Yl.l)  -  . -1(1.1)  ( 3.3 1 

1'6  =  ,T(.\T,2)  —  A{  1.2)  (3.4) 

V6  =  ARCTAX-rrrr  (3.5) 

A  6 

With  this  angle  determined,  the  next  step  was  to  relate  the  horizontal  value  to  a 
value  along  the  zero  fringe  line.  Since  the  horizontal  component  was  the  projection  of 
the  zero  fringe  line  onto  the  horizontal  axis,  the  distance  along  the  zero  fringe  line.  RS. 
was  found  bv: 


cos(  \  6) 

Converting  this  into  centimeters  yields  the  horizontal  component  of  the  fringe  data 

R$(f)  =  RS  x  2  x  Dl  (3.7) 

The  final  step  in  determining  the  fringe  data  was  to  find  the  fringe  value  of  a  vertical 
component.  This  was  done  by  finding  the  contribution  to  the  vertical  component  made 
by  the  zero  fringe  line  and  subtracting  it  from  the  digitized  vertical  component. 

The  vertical  contribution  of  the  zero  fringe  line  .Y8.  was  determined  from: 

T8  =  AT  tan( l'f>)  I  AS) 

This  value  was  subtracted  from  the  vertical  component  of  a  digitized  point  and 
finding  its  projection  along  the  horizontal  axis. 


D}  =*  I  }  1  —  }  S)  cost  t  6) 
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The  fringe  value  was  finally  obtained  by  multiplying  ID?  by  the  ratio  of  fringes  to 
diem/me  units 
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Hence  the  fringe  data  originally  observed  on  the  interferogram  was  converted  into 
a  numeric  form.  RSiJ)  ,  Y5(J>  which  allowed  mathematical  processing. 


E.  SUMMARY 

Digitizing  the  interferogram  was  necessary  to  allow  mathematical  processing  by  the 
computer.  In  the  process,  determining  the  relationship  between  the  fringe  displacement 
and  the  zero  fringe  line  was  required.  Hawns  accomplished  tins,  determining  the  density- 
profile  and  subsequently  the  characteristics  of  the  reduced-density  channel,  using  the 
computer  was  possible.  This  determination  was  conducted  using  two  separate  tech¬ 
niques  which  arc  addressed  in  the  following  chapter. 


IV.  GAUSSIAN  ANALYSIS 

A  Reduced- Density  Channel  was  generated  by  firing  a  C02  laser  beam  into  a  test 
cell  filled  with  a  gas  mixture  of  \2  and  NII3.  Mach-Zehnder  type  interferometers  were 
positioned  at  different  locations  along  the  test  cell  and  the  resulting  interference  patterns 
were  recorded  on  photographic  film.  These  records,  interferogranis,  were  analyzed  using 
Gaussian  Analysis  and  Abel  Inversion  techniques.  This  chapter  is  devoted  to  the  devel¬ 
opment  of  the  Gaussian  technique  and  the  numerical  manipulation  required  to  obtain 
the  characteristics  of  the  reduced-density  channel  from  the  digitized  data  files  using  the 
Gaussian  technique.  The  Abel  technique  is  addressed  in  the  next  chapter. 

A.  THEORY 

The  Gaussian  technique  used  to  determine  the  characteristics  of  a  reduced  density 
channel  was  originally  developed  by  Czuchlewski  for  his  experiment  with  a  test  cell 
containing  N2  slightly  doped  with  SF6  [Ref.  8:  pp.  11-14],  This  technique  will  be  pre¬ 
sented  by  first  considering  a  uniform  square  channel  and  obtaining  a  relationship  be¬ 
tween  the  on-axis  density,  fringe  displacement  and  full  width  at  half  maximum.  The 
square  channel  case  will  then  be  generalized  to  cover  a  non-uniform  channel.  The  result 
will  be  a  relationship  similar  to  the  square  channel  case,  but  applicable  to  the 
interferogranis  obtained  in  the  reduced-density  channel  experiments.  In  this  develop¬ 
ment  a  link  between  the  Gaussian  and  Abel  techniques  will  be  identified  along  with  key 
values  used  in  both  techniques.  These  key  values  are  listed  in  Table  2.  Having  deter¬ 
mined  a  relationship  between  the  on-axis  density,  fringe  displacement  and  full  width  at 
half  minimum  for  a  non-uniform  channel,  its  application  to  the  stored  data  files  will  be 
described.  The  results  of  the  Gaussian  Analysis  for  each  fringe  line  are  presented  in 
Tables  3  through  6. 

1.  Square-Channel 

Consider  a  hypothetical  situation  in  which  one  leg  of  the  interferometer  is 
passing  through  a  square  channel.  This  channel  and  the  resulting  interference  pattern 
are  shown  in  Figure  9.  Designating  the  width  of  the  channel  as  T.  the  optical  path  dif¬ 
ference,  OPD,  between  the  arms  of  the  interferometer  can  be  related  to  the  index  of  re¬ 
fraction  in  the  channel  and  in  the  reference  beam  by: 

r,p [)  =  [>,  _„)  x  r  1 4.  i  j 
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•  nc  is  the  index  of  refraction  in  the  channel. 

•  n,  is  the  index  of  refraction  in  the  reference  arm  of  the  interferometer 

One  fringe  shift  equates  to  a  phase  change  of  2-  lire  fringe 
Figure  9.  AF.  is  the  optical  path  dilferencc  divided  bv  the  wj\  eleug'di.  /.  Sei'M 
in  equation  4.1  yields  a  relationship  between  fringe  shift.  A1  .  and  channel  v. 


AF=  (nc-n0)~ 


The  index  of  refraction  is  related  to  the  mean  polari/ai'ihtv .  ?,  of  a  gas  m  <  (j 
notation  through  the  Lorenz-Lorcntz  relation  (Ref.  14:  p  s"; 
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where 

•  £  is  the  dielectric  constant. 

•  n  is  the  index  of  refraction. 

•  N  is  the  number  of  molecules  per  unit  volume 

•  c  =  n:- 

Making  use  of  the  approximation  that  the  index  of  rcfiuctmn  o!  a  gas  n  near: 
unity,  equation  4.3  can  be  reduced  to: 

n~  —  \ 

Rcarranging  equation  4.4  and  solv  ing  for  the  index  of  refraction  v  :c!dv 

n  =  ( 1  4-  4 -a  A”)  1  i4  : 

Applying  the  binomial  expansion  theorem  to  equation  4  5  and  discarding  tiic  higher  m 
der  terms  yields: 

.>/  =  1  —  2~x.V  4. 


In  order  to  obtain  a  relationship  be*. ‘.seen  the  index  oi' relraet:v:i  and  the  desalts . 
a  connection  between  the  mean  polarizabmr. .  a  .  and  Jcns.ts  must  oe  made  I  :.e  perti¬ 
nent  connection  is  made  through  the  molar  retract:'.  t\ .  A  ' Rot  Id  n  vs1 
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where 

•  R  15  the  ideal  gas  constant 

•  T  ts  the  temperature  in  degrees  Kehm 

•  P  is  the  pressure  of  the  gas 

•  .V,  is  Asagadros  number. 

Sols  mg  equation  4  ”  tor  the  index  ot  retraction,  m  and  using  oiii>  the  first  order 
terms  ot  the  binomial  expansion  theorem  vicids: 


n  »  1  + 


2Rf\ 


Recognizing  that  the  density  of  a  gas.  p  .  can  be  expressed  as: 


and  sols  ing  equation  4.8  for  the  index  of  refraction  in  terms  of  density  yields: 


(4.10) 


From  equations  4.3  and  4.7  one  observes  that  the  mean  polarizability,  y,  and  the  molar 
rcfractivity.  A.  are  constants.  This  is  confirmed  by  Born  and  Wolf  [Ref.  14:  p.  88]  as  they 
demonstrate  that  die  molar  rcfractivity  of  air  varies  less  than  2  percent  between  pres¬ 
sures  of  1  and  42  atmospheres.  In  their  development  of  the  Gaussian  analysis. 
Czuchlewski  et  ah.  'Ref.  8:  p  .11],  regard  the  index  of  refraction  as  a  product  of  the 
density,  p,  and  a  constant.  Hence  equation  4.10  can  be  written  as: 

«  -  1  -  K/i  (4. 1  1 ) 
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Substituting  equation  4.11  into  equation  4.2  yields  a  relationship  between  tire 
fringe  shift  and  the  densitv  of  the  gas  mixture. 


( Krp,  -  K  .pjr 


t\F  = - — - 

/ 

where 

•  pc  is  the  density  of  the  gas  in  the  channel 

•  p„  is  the  density  of  the  gas  at  ambient  conditions 


(4.13) 


Since  it  has  been  determined  that  A.  and  subsequently  K.  is  constant,  equation 


4.13  can  be  written  as: 


A  F  = 


K  nPc-p0) 


(4.14) 


Defining  <5  as  the  ratio  of  the  on-axis  density  to  the  density  at  the  channel  wall, 
allows  equation  4.14  to  be  written  as; 

4„  (<3  -1)KI>„ 

— - - -  (4.1: T 

Hence  equation  4.15  provides  a  relationship  which  allows  one  to  determine  the 
ratio  of  densities  in  a  square  channel  test  cell  by  measuring  the  fringe  shift  and  the  full 
width  at  half  maximum. 


0  =  1- 
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Having  developed  equation  4.16  (or  the  case  of  a  square  channel,  the  next  step 
is  to  generalize  the  procedure  for  a  channel  which  is  indicative  of  the  interlerograms. 

2.  Non-Uniforin  Channel 

Consider  the  interlerogram  sketched  in  figure  10.  This  represents  a  non-uniform 
channel  in  which  the  channel  density  and  index  of  refraction  are  functions  of  position. 
(XA  ).  The  optical  path  difference  at  a  position  along  the  y  axis  can  be  obtained  h\  an 
integral  alone  the  \  axis. 


•Xi 


m 
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OP  [){)-)=  'L_Llnc(X^-noyX 
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where  L  is  large  enough  to  extend  beyond  the  limits  of  the  channel.  Using  the  definition 
of  fringe  shift.  F,  and  leplacing  the  indices  of  refraction  by  there  equivalence  m  terms 
of  densitv  vields: 


W) 


K  Po  rL 

=  ~LJ_Lio(X,\-)-\yX 
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Therefore  the  square  channel  case,  when  expanded  to  cover  a  general  non- 
uniform  channel  yields  an  integral  relationship  between  the  fringe  shift  and  the  ratio  of 
densities. 

The  method  chosen  by  Czuchiewski  to  solve  the  integral  [Ref.  8:  p.  13]  was  to 
pick  a  density  profile  which  was  cylindrically  symmetrical.  The  function  chosen.  S  ( r ). 
was  a  Gaussian  function  and  its  characteristics  were  similar  to  the  actual 
interferoarams. 


Hr)  =  1  -  (1  -  d0)  exp( 


where 


—2.8/- 


(4.19) 


•  <55  is  the  ratio  of  on-axis  density  to  ambient  density. 

•  T  is  the  full  width  of  the  channel  at  the  half  minimum  density. 

In  order  to  solve  equation  4. IS  the  chosen  density  profile  was  converted  from 
cylindrical  to  cartesian  coordinates.  This  conversion  allowed  equation  4.19  to  be  written 


?! 


he  integral  ot  equation  4  .•>  Gaussum  u:;d  :!%  *o.aiu'n  a,;ow  v  t!ic  ounce  'in;: 


to  he  obtained  Iron:: 


M ( })  = 


hi-.'  i  c\: 


Evaluating  equation  4.21  at  Y  =  o  \ieMs  an  expression  lor  the  non-uiuiorm 
channel  density  which  closely  resembles  the  relation  obtained  lor  the  square  channel. 


K^,r  t  , 
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Since  AF  and  T  are  measurable  quantities  on  the  inurferogram.  the  ratio  of 
on-axis  densitv  to  ambient  densitv  can  be  determined  from  the  U  iloume: 


:  _  .  m<))  2.S 

0  r  [  k Po(  •  }  ] 


(4  22) 


Hereafter  the  <5„  will  be  referred  to  as  channel  density. 

B.  EVALUATING  TEST  CELL  PARAMETERS 

In  the  previous  section  of  this  chapter  a  relationship  for  <),  was  obtained  in  terms 
of  the  fringe  shift.  AF  .  full  width  at  half  maximum.  T.  and  some  test  cell  parameters, 
equation  4.23.  Before  determining  the  density  profile,  these  parameters  were  evak.med. 

First  consider  the  parameter  /.  which  represents  the  wavelength  of  the  laser  used  in 
the  interferometer.  This  experiment  used  a  IIE-NE  laser  for  the  interferometer  light 
source  which  has  a  wavelength  of  633  nm  in  vacuum.  This  wavelength  changed  with 
each  gas  mixture  used,  however  it  was  determined  that  the  change,  lnm.  was  not  sig¬ 
nificant.  Hence  the  value  of  /  was  fixed  at  633nm  throughout  these  calculations. 

Next,  consider  the  parameter.  />..  this  represents  the  density  of  the  gas  in  the  refer¬ 
ence  leg  of  the  interferometer.  It  was  assumed  that  the  ambient  conditions  for  each  ex¬ 
periment  were  the  same.  300  K  and  600  torr.  From  the  ideal  gas  law  the  ambient  density 
was  found  to  be: 


/’  60(1/ "60  ,  .0  , 

/>  =  T77-  =  - - —  =  3.1(0  x  10  'nutlcsjan  =  I  .Or  >6  x  |0  mnU\'nic\  < m 
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The  third  parameter,  K.  ,  is  proportional  to  the  molar  refractivity  of  the  gas  mixture. 
(See  equation  4.12).  The  following  describes  the  procedure  used  to  obtain  the  molar 
refractivity  and  subsequently  the  value  of  K.  for  each  gas  mixture. 

Born  and  Wolf  state  that  the  molar  refractivity  of  a  gas  mixture  is  equal  to  the 
weighted  average  of  the  molar  refractivities  of  the  constituents  [Ref.  14:  p.  S9j.  for  the 
NH3-N2  experiment  this  can  be  written  as 


*f\7/3-.V2 


N7/r^  V//3 
XH3 


\2 


(4.24) 


where 

•  -Tv/n-.v:  is  the  molar  refractivity  of  the  gas  mixture 

•  Axm  is  the  molar  refractivity  of  NH3 

•  . ■/ v,  is  the  molar  refractivity  of  N2 

•  \vm  is  the  number  of  NH3  molecules  in  the  test  cell 

•  .Ys,  is  the  number  of  N2  molecules  in  the  test  cell 


Therefore  to  obtain  K  for  a  particular  gas  mixture,  the  molar  refractivity  of  Nil 3 
and  \2  must  be  determined.  The  molar  refractivity,  equation  4.'.  is  a  function  of  the 
index  of  refraction  and  the  density. 

Consider  first  molecular  nitrogen,  N2.  In  order  to  determine  us  molar  refractivity. 
A.  the  index  of  refraction  at  633  nm  was  determined.  This  was  accomplished  by  in¬ 
terpolating  between  the  values  given  for  546.1  nm.  1.O00299S.  and  656.3  nm.  l.ooo2‘>S2. 
at  0'  C  and  76U  torr  [Ref.  15].  The  value  obtained.  1  .<>0029854  was  then  corrected  to 
300K  and  600  torr  using  the  Biot-Arago  law  (Ref.  16:  p.  533). 


">  -  l  P 
1  .0O3o7/  760 


where 

•  n,  is  the  index  of  refraction  of  a  gas  at  temperature,  t  (degrees  Centigradci 

•  n0  is  the  index  of  refraction  of  a  gas  at  0-  C  and  7o0  torr 

•  t  is  the  temperature  of  the  gas  in  degrees  centigrade 

•  P  is  the  ambient  pressure  in  torr. 


Having  performed  the  above 


corrections,  tlie  index  of  refraction  foi  \2  at  "  K 


6i)0  torr  was  determined  to  be  l.noo21444o. 
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The  tabulated  value  of  the  index  of  refraction  lor  NH3  at  the  wavelength  of  sodium 
D  light.  5S9.5  nm.  was  available  from  the  Smithsonian  Tables  [Ref.  16:  p.  533],  Insuifi- 
cient  information  prevented  interpolating  to  obtain  the  value  at  633  nm.  Considering  the 
difference  to  be  quite  small,  the  value  at  5S9.5  nm.  1.000376.  was  used  throughout  these 
calculations.  This  value  was  then  corrected  for  temperature  and  pressure  using  the  same 
procedure  outlined  previously  for  N2.  The  index  of  refraction  at  300 K  and  600  torr  was 
determined  to  be  1.0002700. 

The  next  step  in  determining  the  molar  refractivity  was  straightforward  application 
of  equation  4.7.  This  yielded  a  molar  refractivity  for  N2  as  4.517  cm3  per  mole  and  for 
NH3  5.6SS  cm 3  per  mole. 

Having  determined  the  molar  refractivity  for  N2  and  NH3  equation  4.24  was  used 
to  determine  the  molar  refractivity  of  the  gas  mixtures.  These  values  are  listed  in  fable 
2  along  with  the  values  of  K  .  which  were  determined  from  the  previously  developed  re¬ 
lationship 


(4.11) 


As  will  be  shown  later.  K  is  the  link  between  the  Gaussian  and  Abel  techniques. 

Having  determined  the  numerical  values  of  p  and  K.  the  relationships  for  density, 
equations  4.16  and  4.23,  can  be  reduced  to  the  following: 
or  uniform  square  channel 


Tor  non-uniform  square  channel 
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where  d>  = - - 

/ 

I  able  2  lists  the  value  of  </>  and  other  important  parameters  for  each  gas  mixture 
used  in  the  experiment. 


C.  TEST  CASE 


Thus  far  in  the  development  of  the  Gaussian  analysis  an  algebraic  expression  has 
been  obtained  which  allows  one  to  determine  the  ratio  of  on-a\is  density  to  ambient 
density  for  a  non-uniform  channel,  equation  4.27.  A  test  case  was  generated  to  check 
the  accuracy  of  this  expression  prior  to  applying  it  to  the  actual  intcrferograms. 

The  test  case  was  a  hypothetical  cell  filled  with  600  torr  N2  at  30OK.  Equation  4  21 
was  used  to  develope  a  fringe  pattern  in  which  T  was  equal  to  the  squtue  root  of  2.8  and 
50  chosen  to  be  0.1.  These  conditions  reduced  equation  4.20  to  the  following: 

AF{  Y)  =  5.40  exp(  —  Y2)  (4.28) 

Equation  4.28  was  then  solved  for  the  value  of  Y  which  yielded  one-half  the  maxi¬ 
mum  fringe  shift.  Twice  this  value,  1.67  cm.  corresponds  to  the  full  width  at  half  maxi¬ 
mum.  Evaluating  equation  4.27  at  this  value  of  T.  yielded  a  channel  density  of  0.10.  Thus 
the  equations  developed  provided  self-consistent  results. 

D.  APPLICATION  TO  EXPERIMENTAL  DATA 

Having  determined  that  the  equations  developed  for  a  non-uniform  channel  pro¬ 
vided  self-consistent  results,  each  fringe  line  was  analyzed.  Due  to  the  asymmetry  in  the 
actual  interlerogram.  (See  Figure  11).  three  different  values  of  the  full  width  at  half 
maximum  were  obtained.  These  different  values  correspond  to 

•  T,  the  value  obtained  if  the  right  side  were  the  same  as  the  left  side. 

•  T,  the  value  obtained  if  the  left  side  were  the  same  as  the  right  side. 

•  r3  the  value  obtained  when  the  left  and  right  side  values  are  averaged. 

A  computer  program  was  developed  which  determined  the  value  of  T  and  the  value  of 
AF  and  then  computed  the  value  of  the  channel  density.  The  results  of  this  program  are 
listed  in  Tables  3  through  8. 
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E.  SUMMARY 


The  characteristics  of  the  reduced-density  channel  are  determined  by  the  channel 
density  and  the  full  width  at  half  maximum  at  a  particular  location  along  the  test  cell. 
These  characteristics  have  been  determined  from  a  method  of  Gaussian  analysis  devel¬ 
oped  by  CzuchlcWiki  and  are  listed  in  Tables  3  through  8.  While  this  method  assumes 
the  density  profile  is  Gaussian,  it  provides  a  rapid  means  of  determining  an  approximate 
density.  Recalling  equation  4.27 


» 
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=  1  - 


AF(0) 
1.059 4>V 


(4.27) 


it  is  clear  that  upon  measuring  the  peak  fringe  shift  and  the  full  width  at  half  maximum 
a  rapid  determination  of  the  channel  density  can  be  made. 

In  the  next  chapter  the  Abel  technique  will  be  used  to  obtain  the  same  reduced- 
density  channel  characteristics. 


> 

l 

I 

I 

» 

t 


V.  ABEL  ANALYSIS 


In  the  previous  section  a  method  for  determining  the  ratio  of  on-axis  density  to 
ambient  density  and  width  of  the  reduced  density  channel  was  developed  assuming  the 
density  profile  had  a  Gaussian  shape.  In  this  section  a  method  will  be  developed  which 
makes  no  such  assumption  concerning  the  density  profile,  but  requires  the  profile  to  be 
cylindricallv  symmetric.  This  technique  was  originally  developed  by  Neils  Heinrik  Abel 
in  the  early  1800's  in  his  solution  of  a  mechanics  problem.  In  this  section  the  Abel 
technique  will  be  developed,  applied  to  a  test  case  and  then  applied  to  the  interferometric 
data  obtained  during  the  reduced  density  channel  experiment  involving  NH3-N2  gas 
mixtures.  The  results  of  this  analysis  are  listed  in  Tables  9  through  14  and  provide  a 
determination  of  the  density  profile  which  is  representative  of  the  interferometric  data. 

A.  ABEL  INTEGRAL  TRANSFORMATION 

In  his  short  life  Neils  Heinrik  Abel  made  numerous  contributions  to  the  field  of 
mathmatics  and  science.  The  significance  of  many  of  his  contributions  were  not  recog¬ 
nized  until  more  than  100  years  after  his  death  in  the  late  1820's.  One  area  of  science 
which  has  greatly  benefited  from  his  work  is  the  study  of  plasma  spectroscopy. 

While  attempting  to  solve  a  mechanics  problem  regarding  the  velocity  of  an  object 
along  an  arc  [Ref.  17:  p.  97],  Abel  determined  that  if  a  given  function,  H(a),  can  be  re¬ 
lated  to  another  function  VF  by 


H(a)  = 


o  (a-.t)'' 


then  VF  can  be  obtained  from 


(5.1) 


VF  = 


sin  ,u~ 

71 
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where  0  <,  ii  <  1.  The  full  development  of  this  result  is  available  in  pages  97  through 
103  of  Reference  1 7. 
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As  mentioned  this  integral  transformation  is  widely  used  in  plasma  spectroscopy.  In 
that  application  the  emissivity  of  a  plasma  source  is  determined  from  the  intensity 
measurements.  In  the  development  of  this  emissivity-intensity  relationship  for  use  in 
plasma  spectroscopy.  Orient  [Ref.  18:  p.  176]  found  the  intensity  of  an  optically  thin 
plasma  contained  by  a  cylindrical  column  of  radius  R0  to  be 


r<*o  -  rY‘ 

/(K)  =  2jQ  e(r)dX  = 


c(r)dr 

r  (r2-Y2f 


Using  Abel's  integral  transformation  the  following  determination  of  the  emissivity  was 
obtained  (See  Figure  12) 
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where  — jy—  is  the  derivative  of  I(Y)  with  respect  to  Y.  This  same  method  is  used 
throughout  plasma  physics  in  order  to  determine  the  emissivity  of  the  source  (See  Refs. 
19  and  20  for  examples). 

Typically  the  intensity  is  a  measured  quantity  and  not  a  differentiable  function, 
hence  the  solution  is  accomplished  by  numerical  analysis.  Numerous  methods  of  solving 
equation  5.4  numerically  have  been  identified,  some  of  these  may  be  (bund  in  references 
21  through  25. 

Having  identified  the  Abel  integral  transformation  and  one  of  its  more  widely  used 
forms,  the  next  step  was  to  find  an  expression  between  the  fringe  data  and  the  density 
which  satisfied  the  relationship  expressed  in  equations  5.1  and  5.2  .  Once  identified,  the 
Abel  transformation  was  performed  and  resulted  in  an  expression  lor  the  density  as  an 
integral  function  of  the  fringe  shift. 

B.  APPLYING  THE  ABEL  TRANSFORMATION 

In  the  Gaussian  development  an  integral  expression  between  fringe  shift  and  density 
was  obtained  for  the  case  of  the  non-uniform  channel,  equation  4. IS.  Considering  the 
radial  symmetry  of  the  gas.  defined  by 

r2  =  A'2  +  Y2  (5.5) 


Equation  4,18  was  written  as 
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Substituting  for  JX,  its  equivalent  yielded 
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Performing  the  Abel  Transformation  on  equation  5.7  yielded  an  expression  for  ratio  of 
the  channel  density  to  ambient  density 


8{r)  ~  1  = 
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where  is  the  derivative  of  AF  with  respect  to  Y.  Hence  the  ratio  of  the  densities. 

P(r),  can  be  determined  by  solving  the  right  side  of  equation  5.8.  Attempting  to  solve 
this  analytically  for  any  but  the  simplest  cases  is  extremely  difficult.  Consider  as  an  ex¬ 
ample  the  Gaussian  test  case  of  the  previous  chapter.  Substituting  the  derivative  of  the 
Gaussian  fringe  pattern  into  equation  5.8  produces 
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Thus  even  for  the  situation  where  the  fringe  data  is  given  as  a  differentiable  function, 
numerical  analysis  is  required  to  solve  the  right  side  of  equation  5.8. 

Having  determined  that  the  Abel  Integral  transformation  can  be  applied  to  the  case 
of  the  reduced-density  channels,  the  next  step  is  to  determine  the  method  used  to  obtain 
a  numerical  solution  to  equation  5.8. 

As  mentioned  use  of  the  Abel  integral  transformation  in  determining  the  emissivity 
of  a  plasma  source  is  well  documented.  It  is  believed  that  this  is  one  of  the  first  uses  of 
this  method  in  determining  the  density  profile  from  interferometric  data.  Solution  of  tins 
expression,  like  most  .Abel  integrals  is  best  performed  by  numerical  analysis.  In  this  case 
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an  analytic  solution  is  not  readily  available.  In  the  case  of  plasma  spectroscopy  the  in¬ 
tensity  is  a  measured  quantity,  and  not  a  diffeientiable  function.  The  next  section  de¬ 
scribes  the  method  used  to  determine  the  reduced  density  channel  characteristics  using 
a  numerical  solution  to  the  Abel  integral. 


C.  NUMERICAL  SOLUTION 

A  sketch  of  the  test  cell  is  provided  in  Figure  13.  This  slice  of  the  test  cell  provides 
the  foundation  upon  which  the  Abel  technique  was  developed  for  obtaining  the  reduced 
density  channel  characteristics  of  the  NH3-N2  gas  mixtures. 

In  the  development,  it  was  assumed  that  a  radially  symmetric  density  profile  existed 
after  the  gas  was  energized  by  the  incident  C02  laser  beam.  Further  it  was  assumed  that 
the  region  could  be  divided  into  an  infinite  number  of  concentric  regions  in  which  the 
index  of  refraction  was  constant.  To  illustrate  this  Figure  13  is  di  vided  into  four  regions. 
The  index  of  refraction  at  all  locations  in  region  a  is  the  same.  Likewise  the  index  of 
refraction  for  all  locations  in  region  y  is  the  same  and  so  on  for  regions  >/  and  v.  This  is 
not  to  imply,  however  that  the  index  of  refraction  in  region  a  is  the  same  as  in  regions 
y  ,  )]  or  v.  Having  made  these  assumptions,  an  expression  for  the  optical  path  difference 
along  a  chord.  C„  was  determined  in  terms  of  the  notation  in  Figure  13. 


Q  =  2{(«,  -l)(r2  -  +  («„  -1)[(^ -rf)n  -  (rf  -  rfr"]  +  ....} 

Expanding  this  as  a  series  of  N  regions  yields  the  following 


(5.10) 


C,  =  2  Y  (ftj  -  l){(r2+I  -  rff1  -  ( rj  -  rff'} 


(5.11) 


In  order  to  obtain  an  expression  for  the  index  of  refraction,  the  right  side  of  equation 
5.11  was  manipulated  to  obtain 


Q  -  Y  2 (Hi  -1  ){(Y+,  -  rfy1  -  (rf  -  i'.2)1} 


(5.12) 


In  continuing  with  the  development  of  the  Abel  technique,  a  relation  between  the 
index  of  refraction  and  density  was  required.  This  relationship  was  identified  in  the  de¬ 
velopment  of  the  Gaussian  technique  as  equation  4.11. 


n  —  1  =  Kp 


(4.11) 


Writing  the  previous  equation  in  terms  of  density  and  substituting  the  equivalence  of 
optical  path  difference,  OPD  =  AF  (V)  /,  yielded  the  following 


ATO/  -  £  2 {nj  -  \){{rj+x  ~  r?)Vi  -  {rj  -  rff'} 
2K(r}+x  -  r;2)'/l 


(5.13) 


This  equation  can  be  solved  by  starting  at  the  outer  edge  of  the  test  cell.  X  =  A\  .and 
proceeding  towards  the  center,  X  =  0.  Having  obtained  an  expression  lor  p„  the  next 
step  was  to  determine  the  ratio  of  gas  density  in  the  channel  to  the  ambient  density.  This 
ratio  was  obtained  by  scaling  the  density  in  the  channel.  Recall  that  the  density  along 
the  axis  was  determined  to  be  less  than  the  density  along  the  walls  of  the  test  cell.  De¬ 
fining  a  scale  comparable  to  that  of  the  Gaussian  technique.  was  set  equal  to: 

4„=l-~  (5.14) 

One  can  observe  that  from  equation  5.14,  the  value  of  <)„  is  1  at  the  wall  of  the  cylinder. 
pc<$p0  ,  and  decreases  at  positions  close  to  the  center  of  the  channel. 

Having  determined  a  relationship  for  the  density  in  the  channel,  a  computer  pro¬ 
gram  was  developed  which  extracted  the  density  profile  of  the  channel  from  the  fringe 
data.  This  method  was  then  checked  using  the  Gaussian  test  case  developed  for  the 
Gaussian  technique.  Application  of  this  test  case  to  the  Abel  technique  is  described  in 
the  next  section. 

D.  TEST  CASE 

A  method  of  obtaining  the  density  profile  of  the  test  cell  using  the  Abel  technique 
was  developed  in  the  previous  section.  This  method  was  checked  using  the  fringe  pattern 
described  by  equation  4.2S  and  found  to  provide  reults  ol  d,  which  agreed  with  test  re¬ 
sults  to  within  two  decimal  places;  the  results  off  agreed  to  within  o.07  cm.  Table  15 
lists  the  Abel  and  test  results  lor  <),  and  F  for  each  gas  mixture  used.  As  an  example  of 


the  resultant  density  profile,  the  Abel  determined  profile  for  <><><»  torr  \2  obtained  from 
the  fringe  pattern  described  by  equation  4.2S  is  plotted  in  figure  14.  1  lie  value'-  used 
to  generate  the  plot  are  provided  in  Table  1<>.  I  lasing  determined  the  Abel  method 
produced  results  winch  agreed  with  the  test  cases,  the  density  profile  for  each  limge  line 
digitized  was  determined. 

E.  APPLICATION  TO  EXPERIMENTAL  DATA 

After  determining  that  the  Abel  technique  produced  satisfactory  results,  each  fringe 
line  was  analyzed  to  determine  the  density  profile.  In  doing  so  it  became  clear  that  the 
slanted  interferogram.  characteristic  of  the  M13-N2  gas  mixtures,  presented  a  problem 
that  was  not  encountered  with  the  test  case. 

In  the  test  case,  the  fringe  lines  were  divided  into  two  parts  and  the  density  profile 
of  each  part  matched  at  the  origin.  In  the  case  of  the  actual  interferograms  the  difference 
in  fringe  patterns  on  either  side  of  the  division  line  precluded  a  match  at  the  origin.  This 
difference  is  consistent  with  the  development  of  the  Abel  technique,  however  it  violates 
the  assumption  that  the  density  profiles  were  cylindrically  symmetric.  At  this  point  the 
assumption  was  made  that  each  part  of  the  fringe  line  represented  half  of  a  symmetric 
fringe  line.  The  value  of  <)„  to  be  used  in  the  analysis  of  the  results  was  an  average  of 
the  values  obtained  for  both  parts  of  the  fringe  line.  The  same  procedure  was  used  for 
determining  the  value  of  T.  The  results  of  Abel  analysis  on  each  fringe  line  are  listed  in 
Tables  9  through  14.  The  values  of  <5,  are  provided  for  both  sides  of  the  division  line  as 
well  as  the  average  between  the  two  sides.  These  three  values  can  then  be  thought  as 
defining  the  minimum,  maximum  and  most  likely  value  of  the  on-axis  density  and 
channel  width. 

F.  SUMMARY 

The  Abel  technique  provided  a  means  to  determine  the  density  profile  of  a  gas 
mixture  heated  by  a  C02  laser  beam.  From  this  density  profile  a  determination  of  the 
reduced  density  channel  characteristics  was  possible.  In  the  next  chapter  results  of  the 
the  Gaussian  and  Abel  techniques  will  be  used  to  determine  relationships  between  the 
channel  characteristics  and  the  variables  of  the  experiment,  energy  and  NII3  pressure. 


VI.  RESULTS 


Interferograms  made  shortly  alter  a  CG2  laser  energized  a  test  cell  containing  vari¬ 
ous  mixtures  of  N2  and  XH3  have  been  analyzed  using  two  different  methods.  These 
methods,  Gaussian  and  Abel,  have  been  described  in  detail  in  the  preceding  chapters.  In 
this  section  a  comparison  of  the  two  methods  and  some  general  observations  will  be 
made  concerning  channel  formation. 

A.  GAUSSIAN- ABEL  COMPARISON 

A  summary  of  the  experimental  results  is  provided  in  Tables  17-22.  Included  in 
these  tables  is  a  comparison  between  the  Gaussian  and  Abel  determined  values.  The 
closest  agreement  between  the  two  methods  can  be  found  in  Table  19.  In  that  Table  the 
results  of  the  experiment  for  a  test  cell  containing  a  gas  mixture  of  15  torr  NI13  are 
listed.  At  112  cm  along  the  test  cell  and  for  a  laser  energy  of  84  joules,  the  Gaussian  and 
Abel  determined  values  of  <)„  differ  by  1.7  percent.  Under  the  same  conditions  the  values 
determined  for  T  differ  by  3.6  percent.  This  fringe  pattern  and  the  Gaussian  curve  ob¬ 
tained  by  using  values  of  AF  and  T  obtained  from  the  fringe  data  are  plotted  in  Figure 
15.  As  expected  the  two  methods  yield  similar  results  when  the  fringe  pattern  has  a 
Gaussian  shape. 

In  some  cases,  however,  the  fringe  data  are  not  Gaussian  shaped  and  the  analytic 
analysis  is  not  expected  to  be  very  accurate.  Consider  the  results  of  5  torr  NII3  listed 
in  Table  17.  At  51  cm  and  for  an  energy  of  31  joules,  the  values  of  and  T  obtained 
by  the  Gaussian  and  Abel  methods  differ  appreciably.  For  b,  the  two  analysis  methods 
yield  results  which  differ  by  21.1  percent;  the  salues  for  T  differ  by  nearly  331  percent! 
This  fringe  line  and  the  corresponding  Gaussian  curve  are  plotted  in  Figure  16.  It  is 
apparent  from  the  differences  between  the  curves  that  the  fringe  pattern  is  not  Gaussian 
shaped.  Recalling  that  the  Gaussian  method  uses  the  maximum  fringe  shift  and  the  full 
width  at  half  maximum  to  determine  the  value  of  and  that  the  Abel  method  uses  all 
the  points  digitized  along  the  fringe  line,  it  is  easily  understood  why  the  two  methods 
yield  different  results  for  non-Gaussian  shaped  fringe  patterns. 

In  general  it  can  be  stated  that  tiie  smaller  the  difference  between  the  Gaussian  and 
Abel  determined  values  of  5,  and  T  the  more  nearly  the  fringe  line  approximates  a 
Gaussian  shape.  A  review  of  Tables  17-22  reveals  that  the  Gaussian  determined  T  and 
the  Abel  determined  f  differ  by  more  than  2o  percent  in  all  but  8  of  the  fringe  lines. 


The  corresponding  values  of  <i,3  in  these  instances  differed  by  less  than  5.3  percent. 
Hence  less  than  20  percent  of  the  experimental  data  was  Gaussian  shaped.  For  this 
reason  analysis  of  the  reduced-density  channel  formation  will  be  based  on  the  results 
obtained  from  Abel  method. 

The  Gaussian  method  was  developed  to  provide  a  quick  assessment  of  the  channel 
formation.  From  these  experimental  results  it  can  be  stated  that  for  fringe  patterns  of 
Gaussian  shape,  it  yielded  results  very  similar  to  those  obtained  from  the  Abel  method. 
For  non-Gaussian  shaped  fringe  patterns,  it  can  be  stated  that  it  yields  values  of  5, 
within  31  percent  of  the  Abel  method.  As  for  the  values  of  F  obtained,  the  Gaussian 
method  provided  results  which  were  always  greater  than  or  equal  to  the  Abel  results. 
Hence,  for  this  experiment  the  Gaussian  method  identifies  the  largest  channel  width 
possible. 

B.  ACCURACY  OF  RESULTS 

In  order  to  provide  a  determination  of  the  accuracy  of  the  entire  process  from  pho¬ 
tographic  enlargement  through  the  determination  of  50,  a  series  of  error  estimates  were 
made.  In  these  calculations  it  was  assumed  that  the  largest  probable  error  in  digitizing 
a  fringe  line  would  be  10  percent  of  the  actual  value.  Using  the  fringe  pattern  described 
by  equation  4.2S  two  new  fringe  patterns  were  generated  each  dilfering  from  the  original 
by  10  percent.  One  of  the  fringe  patterns  generated  consistently  differed  by  lu  percent 
more  than  the  original  pattern:  the  second  pattern  generated  consistently  differed  by  10 
percent  less  than  the  original  pattern.  The  results  of  the  Abel  determined  values  of  b  at 
each  gas  mixture  are  plotted  in  Figure  17.  As  can  be  observed  a  difference  of  ±  10 
percent  in  the  fringe  pattern  results  in  a  ±  0.09  differnce  in  <>„.  Since  the  assumption  of 
10  percent  error  incorporates  the  error  introduced  in  the  photographic  enlargement 
process  and  covers  erratic  digitizing  errors,  there  is  confidence  that  the  values  of  b  are 
accurate  to  within  ±  0.09. 

C.  ANALYSIS  OF  RESULTS 

The  Abel  method  is  more  dependent  on  the  shape  of  the  entire  fringe  line  than  the 
Gaussian  method.  To  gain  insight  into  the  meaning  of  the  results  listed  in  fables  l~-22. 
several  plots  of  the  Abel  determined  values  of  a  and  T  were  made.  The  results  could 
be  plotted  in  numerous  ways,  for  example  <>.  as  a  function  of  MI3  concenttaiion.  laser 
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energy,  distance  from  the  entrance  of  the  cell.  etc.  The  methods  chosen  to  represent  the 
results  are  60  as  a  function  of  laser  energy.  Figures  18-22.  and  0,  as  function  of  distance 
from  the  test  cell  entrance.  Figures  23  and  24.  These  two  methods  will  provide  a  means 
to  compare  the  results  to  the  results  of  the  Sl'o  reduced-density  channel  experiment. 

1.  Channel  Density 

The  following  observations  were  made  from  the  plots  of  o  as  a  function  of  en¬ 
ergy  and  with  NH3  concentration  and  distance  from  cell  entrance  constant  parameters 
(Figures  IS  through  22). 

At  51  cm  from  the  entrance  of  the  test  cell  {  See  Figures  IS  and  19): 

•  The  plots  for  5  torr  and  200  torr  NII3  indicate  that  there  is  a  slight  decrease  in  5. 
as  energy  is  increased.  One  could  also  fit  the  data  with  a  constant  value  lor  0.  of 
0.73  +  *0.05  for  5  torr  and  0.S3  ±  0.03  for  200  torr.  Hence  an  equally  valid  inter¬ 
pretation  of  the  results  is  that  there  is  no  apparent  energy  dependence  on  the 
channel  depth  (  Figures  IS  and  19). 

•  For  15  torr  (See  Figure  IS)  a  pronounced  dependence  of  a.  on  the  energy  is  ob¬ 
served.  channel  density  decreases  with  increasing  energy.  At  low  energies  o.  is  ob¬ 
served  at  values  around  0.60.  while  at  higher  energies  values  around  O.?o  are 
observed. 

•  For  50  torr  the  channel  density  decreases  with  increasing  energy  through  the  range 
of  0  to  43  joules.  At  43  joules  there  appears  to  be  a  minimum  \aluc  and  then  the 
channel  density  begins  to  increase.  This  minimum  will  later  be  sompated  to  the 
minimum  predicted  by  the  computer  model  for  the  S I  o  experiment 

At  112  cm  from  the  entrance  to  the  test  cell  (  See  Figures  2<*  and  21 ): 

•  For  5  torr  there  appears  to  be  little  effect  on  n.  when  energy  is  increased.  1  iie  re¬ 
sults  could  be  fitted  with  a  constant  6.  of  0.“3  ±  <M»7  ( [  igure  2<)i. 

•  For  15  torr  a  distinguishable  dependence  of  i)„  is  observed,  channel  densitx  de¬ 
creases  as  energy  increases  (Figure  2t>). 

•  For  50  torr  one  can  observe  remnants  of  the  minimum  identified  at  51  cm  (Ficure 
21). 

At  221  cm  from  the  entrance  to  the  test  cell  (See  f  igure  22): 

•  For  5  and  15  torr  NH3  a  slight  decrease  in  channel  densitx  is  observed  with  in¬ 
creasing  energy. 

•  For  50  and  200  torr  NI13  all  the  laser  energy  is  absorbed  at  the  beginning  of  the 
cell  and  thus,  there  are  no  channels  at  221  cm. 

The  results  obtained  for  0,  are  consistent  and  viable.  1  no  following  identify  treats 
supporting  this 

•  Plots  of  ,5  versus  energy  (Figures  IS  through  22)  ;. icld  smooth  uir  es;  there  are  no 
discontinuities. 
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•  The  values  obtained  for  J,  under  similar  conditions  yield  similar  results.  An  exam¬ 
ple  is  the  results  obtained  for  15  torr  at  51  cm  and  S4  joules  ( Figure  IS).  The  values 
of  <)„  are  close  and  are  within  the  statistical  error  bar. 

To  provide  further  comparisons,  representative  data  points  were  selected  and 
plotted  as  a  function  of  position  along  the  test  cell  in  Figures  2'  and  24.  The  energies 
were  selected  to  provide  a  means  of  comparison  between  different  concentrations  of 
N‘H3.  Among  the  more  prominent  features  of  each  plot  in  Figure  23  are 

•  As  the  distance  from  the  entrance  of  the  test  cell  is  increased  the  value  of  <3.  in¬ 
creases.  This  finding  agrees  intuitively  with  the  idea  that  the  laser  beam  is  contin¬ 
uously  absorbed  as  it  travels  further  from  the  entrance.  Hence,  there  is  less  energy 
available  at  the  end  of  the  channel  than  at  the  beginning  and  we  expect  weaker 
channels  at  the  end  of  the  cell.  The  situation  for  5  torr,  45  joules  and  at  distances 
less  than  130  cm  from  the  entrance  of  the  cell  may  be  more  complicated. 

•  As  the  concentration  of  X II 3  increases,  the  length  ol  the  channel  decreases.  This 
is  evidenced  by  the  fact  that  no  channels  exist  in  2<»/  torr  at  1 12  cm  and  none  exi't 
at  221  cm  for  50  torr.  It  also  agrees  intuitively  with  the  idea  that  as  the  number 
of  absorbers  per  unit  volume  increases,  the  distance  a  beam  travels  before  dissi¬ 
pating  decreases. 

•  As  the  concentration  of  NH3  increases  from  5  to  15  torr.  channel  density 
decreases. This  is  true  for  all  energies  plotted  and  is  consistent  with  the  findings  of 
the  N2-SF6  experiment.  In  short,  the  higher  the  concentration,  the  deeper  the 
density  channel  (Figure  2<j). 

•  For  a  laser  energy  of  45  joules  the  channel  produced  in  5  torr  Nil.'  is  relative!;* 
constant  over  the  range  of  5<>  to  221  cm  with  a  value  of  <3  of  o.“h  -  <u>5.  1  h;s  is 
also  true  for  75  joules,  for  which  0.  remains  0.6$  r  ‘MU.  This  is  similar  to  the  re¬ 
sults  obtained  for  the  lowest  concentration  of  SFo  |Ref.  S:  p.  2“j.  for  which  rela¬ 
tively  long  shallow  channels  were  observed. 

2.  Channel  Width 

The  channel  widths  for  each  concentration  at  selected  energies  arc  plotted  in 
Figure  24.  The  selected  energies  arc  the  same  used  in  Figure  23  and  allow  a  a  comparison 
between  values  at  the  same  energy,  Unlike  the  plots  of  o,  versus  distances,  very  little 
consistent  information  can  be  garnered  from  Figure  24,  One  identifiable  ohuiuoteristic 
can  be  observed  by  considering  the  5  and  15  torr  concentrations  at  51  cm  and  1 12  cm. 
In  these  instances  the  channel  width  increases  as  energy  increases.  The  opposite  appears 
true  for  20<>  torr  NII3.  F  decreases  as  energy  increases. 
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D.  COMPARISON  WITH  SF6 

At  present  there  is  no  quantitative  method  K. r  comparing  the  characteristics  of  the 
reduced-density  channel  obtained  in  the  N1I3-N2  experiment.  For  that  reason  the  results 
will  be  compared  with  the  results  of  the  N2-SF6  reduced-density  channel  experiment 
[Ref.  8:  p.  27[. 

The  results  of  the  N2-SF6  and  the  N2-NH3  experiments  were  very  similar.  The 
major  differences  being  that  NII3  was  substituted  for  SF6  in  the  tc^t  cell  as  the  !»'.<> 
ixm  radiation  absorber  and  the  measurements  were  taken  at  dtlTerem  distances.  A  com¬ 
puter  model,  CHANL,  developed  by  Naval  Research  Laboratory  scientists  was  used  to 
predict  the  results  of  the  N2-SF6  experiment  [Ref.  8:  p.  S|.  A  signlicant  product  of  that 
model  is  plotted  in  Figure  25.  This  figure  shows  the  predicted  channel  density  lor  a  laser 
energy  of  50  joules  as  a  function  of  distance  from  the  entrance  of  the  test  cell.  The 
N2-SF6  experimental  results,  shown  in  Figure  26.  correlate  '.veil  with  the  results  of  the 
model.  A  trend  of  increasing  absorber  concentration  resulting  in  shorter,  deeper  channels 
is  apparent  in  both  plots.  When  comparing  Figures  25  and  26  with  the  N2-N11.3  results 
shown  in  Figure  23,  the  same  pattern  clearly  exists  for  5,  15  and  50  torr  NI13.  Less 
obvious  is  the  correlation  with  200  torr  NH3.  While  only  two  data  points  are  available, 
one  can  assume  that  the  laser  energy  has  been  almost  totally  absorbed  before  it  reached 
51  cm.  This  assumption  is  supported  by  Figures  26  and  27. 

Figure  27  is  another  product  of  the  CHANL  code.  This  plot  shows  the  ellect  on  6. 
at  0.4  torr  SF6  as  energy  is  increased.  Significant  is  the  minimum  0  which  occurs  at  5<> 
joules  and  the  sharp  change  in  slope  at  165  cm.  Evidence  of  these  changes  in  slope  are 
evident  in  both  the  N2-SF6  (Figures  28  and  29)  and  N2-NH3  (Figures  1)  and  21)  ex¬ 
perimental  results. 

Based  on  the  above  observations,  it  can  be  concluded  that  the  N2-NII3  results  ex¬ 
hibit  the  same  qualitative  behavior  as  the  N2-SF6  results  and  the  CHANL  code  pred¬ 
ictions. 
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E.  CONCLUSIONS 

The  C02  laser  beam  created  reduced-density  channels  in  test  cells  containing  a  gas 
mixture  of  N2  and  NT  13.  These  channels  were  analyzed  and  found  to  be  dependent  on 
many  parameters  including  the  concentration  of  NH3  and  the  laser  energy.  As  expected 
it  was  found  that  the  depth  of  the  channel  decreased  as  the  the  beam  travelled  along  the 
test  cell.  It  was  also  determined  that  for  a  fixed  laser  energy,  the  depth  of  the  channel 
and  its  length  varied  inversely  with  the  concentration  of  NT  13.  That  is.  shorter,  and 
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deeper  channels  were  observed  ac  higher  concentrations  of  N1I3:  long,  shallow  channels 
were  observed  at  lower  concentrations.  This  qualitatively  agrees  with  the  experimental 
results  and  computer  predictions  for  the  N2-S16  experiment  [Ref.  8:  p.  2~\. 

b  ran  also  he  concluded  that  the  Gaussian  metnod  of  analysis  provides  a  quick 
means  of  determining  the  characteristics  of  a  reduced-density  channel.  1  he  accuracy  cf 
the  method  is  dependent  upon  the  shape  of  the  fringe  pattern.  For  Gaussian  shaped 
patterns  the  results  are  highly  accurate.  For  non-Gaussian  shapes  values  for  o,  differed 
from  the  more  accurate  Abel  method  by  as  much  as  31  percent.  The  channel  width  de¬ 
termined  differed  by  as  much  as  330  percent.  Despite  these  large  deferences  the  Gaussian 
method  provided  a  quick  means  of  determining  to  within  .  1 5  and  finding  the  maxi¬ 
mum  value  of  channel  width. 

To  summarize,  the  analysis  of  reduced-density  channels  in  gas  mixtures  of  N2  and 
NH3  support  the  N2  and  SF6  experimental  results.  The  channels  obtained  in  both  ex¬ 
periments  should  be  suitable  for  studying  the  propagation  of  intense  electron  beams 
through  the  atmosphere,  and  thus  provides  a  foundation  for  determining  the  role  these 
beams  will  have  in  the  Strategic  Defense  Initiative. 
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APPENDIX  A.  TABLES 


Table  1.  COORDINATES  OF  ACTUAL  DATA  AND  ENLARGED  DATA 


Table  2.  KEY  VALUES  USED  IN  ABEL  AND  GAUSSIAN  TECHNIQUES 


Pressure  (torr) 

Density  (x  E18  molecules 
per  cubic  cm) 

Molar 
Refractivity 
(cubic  cm 
per  mole) 

K 

(x  C-23 
cubic 
cm  per 
mole¬ 
cule) 

<p(c7/r‘j 

NH3 

N2 

NH3 

N2 

0 

600 

0.00 

19.06 

4.517 

1.1247 

3.3S7 

5 

595 

.1583 

18.84 

4.526 

3.394 

15 

5S5 

4.75 

18.53 

4.546 

1.132 

3.409 

5o 

550 

15.83 

17.42 

4.614 

1.149 

3.461 

200 

4oo 

6.33 

12.67 

4.0(H> 

deb 

3.6"9 

Table  5.  GAUSSIAN  RESULTS  FOR  15  TORR  NH3/585  TORR  N2 


Table  6.  GAUSSIAN  RESULTS  FOR  50  TORR  NH3/550  TORR  N2 


DA  r.\ 
FILE 

FRINGE 

D1ST 

FROM 

SOURCE 

(cm) 

ENERGY 

(joules) 

GAUSSIAN  ANALYSIS 
(ON-AXIS  DENSITY) 

Full 

Left 

Ave 

Test 

I 

NA 

NA 

.124 

.124 

.124 

.124 

S77B5 

5 

112 

63 

.816 

.847 

.768 

.810 

S77B4 

4 

112 

H 

.814 

.815 

.813 

.814 

S77A6 

6 

51 

63 

.4S4 

.513 

.452 

.483 

877A5 

5 

51 

63 

Em 

.490 

.363 

.429 

S7SB7 

7 

112 

76 

.761 

.790 

.724 

.758 

878  B  8 

8 

112 

76 

.769 

.795 

.737 

.767 

S7SA7 

*T 

/ 

51 

76 

.47S 

.432 

.518 

.476 

878A6 

6 

51 

76 

.515 

.477 

.547 

.513 

879B6 

6 

112 

72 

.797 

.777 

.813 

.796 

879B5 

5 

112 

72 

.761 

.784 

.759 

S79A5 

5 

51 

72 

.521 

.521 

.521 

.521 

879A4 

4 

51 

72 

.554 

.570 

.537 

.554 

SS0B7 

7 

112 

73 

.785 

.770 

.79S 

.784 

880B6 

6 

112 

73 

.765 

.752 

.777 

.765 

Table  7.  GAUSSIAN  RESULTS  FOR  50  TORR  NH3/550  TORR  N2 


DATA 

FILE 

FRINGE 

DIST 

FROM 

SOURCE 

(cm) 

ENERGY 

(joules) 

GAUSSIAN  ANALYSIS 
(OX-AXIS  DENSITY) 

Full 

Left 

Right 

Ave 

Test 

1 

XA 

XA 

.124 

.124 

.124 

.124 

973B2 

■) 

112 

12 

.S30 

.S56 

.793 

.826 

973B1 

1 

112 

12 

.854 

.898 

.740 

.830 

973A4 

4 

51 

12 

.641 

.610 

.666 

im— 

973A5 

5 

51 

12 

.683 

.665 

.700 

.683 

974B5 

5 

112 

39 

mm 

.8S2 

.713 

.772 

974B2 

2 

112 

39 

.718 

mm 

.653 

.71  1 

974A7 

7 

51 

39 

.439 

.438 

.441 

.439 

9  74 A  6 

6 

51 

39 

.417 

.414 

.419 

.417 

976B3 

112 

42 

.720 

.774 

.632 

.701 

976B2 

2 

112 

39 

.703 

.778 

.557 

.661 

976A5 

5 

51 

42 

.439 

.352 

.506 

.422 

976A4 

4 

51 

42 

.345 

.259 

.413 

n 

977B5 

5 

112 

43 

.7S0 

.830 

.688 

.751 

977B4 

4 

112 

43 

.739 

.797 

.635 

.716 

3 

51 

43 

.452 

.727 

Em 

977A2 

2 

51 

43 

.567 

.400 

.661 

.530 

Table  8.  GAUSSIAN  RESULTS  FOR  200  TORR  NH3/400  TORR  N2 
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Table  9.  ABEL  RESULTS  FOR  5  TORR  NH3/595  TORR  N2 


DATA 

FILE 

FRINGE 

DIST 

FROM 

SOURCE 

(cm) 

ENERGY 

(joules) 

ABEL  ANALYSIS 
(ON-AXIS  DENSITY) 

Left 

Right 

A  ve 

Test 

1 

NA 

NA 

.102 

.102 

.102 

966A4 

4 

51 

31 

.662 

.733 

.697 

9 66 A  5 

5 

51 

31 

.638 

.773 

.705 

966  B4 

4 

112 

31 

.856 

.752 

.804 

966B5 

5 

112 

31 

.889 

.724 

.806 

964A7 

7 

51 

75 

.713 

.621 

964 A  6 

6 

51 

75 

.716 

.609 

.663 

mm 

J 

112 

75 

.759 

.581 

.670 

964  B  2 

2 

112 

75 

.799 

.565 

.622 

9640 

1 

221 

75 

.5S5 

.840 

.713 

963A4 

4 

51 

45 

.741 

.716 

.728 

963A5 

5 

51 

45 

.S40 

.764 

.802 

963  B  2 

2 

112 

45 

.829 

.608 

.719 

f 

3 

112 

45 

.806 

.(>39 

.723 

963  C  2 

2 

221 

45 

.787 

.894 

.840 

962A4 

4 

51 

57 

.676 

.('SO 

.('78 

962A5 

5 

51 

57 

.731 

.591 

.66 1 

2 

112 

57 

.828 

.774 

.801 

962  B  5 

5 

112 

57 

.864 

.647 

.756 

9620 

1 

221 

57 

.S74 

.860 

.867 

Table  10.  ABEL  RESULTS  FOR  15  TORR  NH3/585  TORR  N2 


Table  14.  ABEL  RESULTS  FOR  200  TORR  NH3/400  TORR  N2 


DATA 

FILE 

FRINGE 

DIST 

FROM 

SOURCE 

(cm) 

ENERGY 

(joules) 

ABEL  ANALYSIS 
i. ON-AXIS  DENSITY) 

Left 

Right 

Ave 

Test 

1 

NA 

NA 

.179 

.179 

.179 

S74A6 

6 

51 

1  IS 

.8  SO 

.833 

.857 

8  74 A  5 

5 

51 

1  IS 

.850 

.848 

.849 

S76A5 

5 

51 

203 

.820 

.SO0 

.810 

876A4 

4 

51 

203 

.700 

.855 

.777 

Table  15.  SUMMARY  OF  TEST  CASE  RESULTS 


Pressure  (torr) 

K  (xE-23 
cubic  cm 
per  mole¬ 
cule) 

Channel  Density  (o0  ) 

Full  Width  at  Half  Min; 
f  (cm) 

NH3 

N2 

Abel 

Te^t 

Abel 

Test 

0 

600 

1.1247 

.102 

.100 

1.60 

1 .67 

5 

mam 

1.127 

.102 

.105 

1.60 

1.67 

15 

5S5 

1.132 

.106 

.109 

1.60 

50 

550 

1.149 

.121 

.124 

1.60 

1.67 

200 

400 

1.222 

.179 

.182 

1.60 

Table  16.  DENSITY  PROFILE  OF  600  TORR  N2  (ABEL) 


Left  Side 

Right  Side 

Dist  from  Origin 
(cm) 

On-Axis  Density 

Dist  from  Origin 
( cm ) 

2.20 

1.00 

2.20 

1.56 

0.95 

1.56 

1.30 

O.SS 

1.30 

1.14 

0.79 

1.14 

1.01 

0.70 

1.01 

0.S9 

0.62 

0.89 

0.7S 

0.54 

0.78 

0.66 

0.47 

0.66 

0.55 

0.37 

0.55 

0.43 

0.2S 

0.43 

0.26 

0.22 

0.26 

0.00 

0.10 

0.00 

4<> 


Table  17.  SUMMARY  OF  RESULTS  FOR  5  TORR  NH3 


Energy 

(joules) 


Dist 

from 

Source 

(cm) 


ON-AXIS  DENSITY  io  ) 


Full  Width  at  I  lall' Mm  ( I~ 


Abel 

± 

Gauss 

± 

Abel 

± 

Gauss 

4- 

.71 

.06 

.86 

0.0 

21.1 

.36 

.03 

1.55 

.06 

.77 

.05 

.82 

.02 

5.2 

.73 

.17 

1.39 

.13 

.68 

.06 

.SO 

.03 

17.6 

.50 

.02 

1.48 

.26 

.65 

.09 

.68 

.02 

4.6 

.79 

.05 

1.11 

.04 

.SI 

.08 

.85 

0.0 

4.9 

.48 

.07 

.SI 

.03 

.73 

.11 

.79 

.03 

S.2 

.SO 

.05 

1.17 

.24 

.78 

.09 

.82 

.02 

5.1 

.88 

.17 

1.41 

.25 

.69 

.12 

.74 

.01 

7.2 

.93 

.11 

1.19 

.<  M 

.85 

.05 

.SS 

.04 

3.4 

.42 

NA 

.78 

NA 

.87 

.07 

.90 

.02 

3.4 

.45 

.03 

.96 

.18 

.72 

.12 

.82 

.05 

13.9 

.45 

NA 

.87 

NA 

33<).6 

90.4 

196.0 

4Q.5 

68.8 


i.3 


60.2 

28.0 

85.7 

46.8 


Table  IS.  SUMMARY  OF  RESULTS  FOR  15  TORR  NH3 


Energy 

(joules) 


Dist 

from 

Source 

(cm) 


ON-AXIS  DENSITY  (d  ) 


Full  Width  at  Half  Min  (T  ) 


Abel 

±  Gauss 

± 

Abel 

+ 

Gauss 

+ 

.58 

.06  .62 

.01 

6.9 

.77 

.05 

1.14 

.03 

.63 

.11  .64 

.02 

1.6 

.79 

.02 

1.04 

.07 

.36 

.11  .36 

.01 

O.o 

.83 

.01 

.95 

.04 

.28 

.18  .34 

.09 

21.4 

.73 

.08 

.89 

,i  i5 

.32 

.(>S  .36 

.07 

12.5 

.95 

.07 

1.21 

.16 

.24 

.16  .24 

.03 

o.O 

.88 

.03 

.94 

0.0 

1 V  v  vv \ 


pjr7\XTk*  ",mjr*.r^ir”  \r* 


Table  19. 


SUMMARY  OF  RESULTS  FOR  15  TORR  NH3 


Dist  ON- AXIS  DENSITY  (S  )  Full  Width  at  Half  Min  (T  ) 

from 

Source 

(cm) 


Abel  |  ±  Gauss  ±  A%  Abel  +  Gauss 


1.17 


.91 


1.15 


NA  59.7 


NA 


ON-AXIS  DENSITY  (o  ) 

Abel 

± 

Gauss 

± 

eh 

.66 

.04 

.72 

.02 

9.1 

.66 

.08 

.69 

.01 

4.6 

.62 

.09 

.65 

0.0 

tm 

.55 

.04 

.64 

.01 

16.4 

.58 

.04 

0.0 

1.7 

.90 

.09 

.91 

NA 

m 

.85 

.01 

.89 

NA 

4.7 

.77 

.05 

.84 

NA 

EB 

Table  20.  SUMMARY  OF  RESULTS  FOR  50  TORR  NH3 


Dist 

from 

Source 

(cm) 


ON-AXIS  DENSITY  (6  ) 


+ 

Gauss 

-f 

.06 

.66 

.03 

.07 

.43 

.02 

.18 

.39 

.07 

.10 

.60 

.05 

.10 

.46 

.04 

■a 

.54 

.02 

.06 

.61 

NA 

.06 

,5o 

.03 

Full  Width  at  Half  Min  (T  ) 


Abel  + 


Gauss 


.99 


1.15 


1.49 


1.16 


.94 


,9S 


Reduced  Density  Channel  Experiment 
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(wo)A 


FRINGE  SHIFT 


DISTANCE  (cm) 


Figure  13.  Cylindrical!}'  Symmetric  Test  Cell  Containing  NH3-N2  Gas  Mixture 
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Figure  19.  Variation  of  Channel  Density  with  Laser  Energy  at  51  cm  from  Entrance 
of  Test  Ceil  for  50  torr  NH3  (bottom),  200  torr  (top). 
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Figure  26.  Variation  of  Channel  Density  with  Distance  from  the  Test  Ceil  Entrance 
for  0.2.  0.4.  0.8  and  2.0  torr  ST6. 
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Director.  Defense  Advanced  Research 
Project  Aacncv 

ATTN:  LCOL  R.A.  Guilickson 
1400  Wilson  Blvd. 

Arlington.  VA  22209-2308 

Defense  Advanced  Research  Project  Asencv 
ATTN:  Maj.  G.P.  Lasche 
1400  Wiison  Blvd. 

Arlington.  VA  22209-2308 

Defense  Advanced  Research  Project  Asencv 
ATTN:  Dr.  Shen  Shey 
Directed  Energy  Office 
1400  Wilson  Blvd. 

Arlington.  VA  22209-230S 

Los  Alamos  National  Laboratory 

ATTN:  Dr.  Randolph  Carlson 

Mail  Stop  P942 

P.O.  Box  1663 

Los  Alamos,  NM  87545 

Los  Alamos  National  Laboratory' 

ATTN:  Dr.  David  Moir 
Mail  Stop  P942 
P.O.  Box  1663 
Los  Alamos,  NM  87545 

Los  Alamos  National  Laboratorv 

ATTN:  Dr.  T.  Neal 

Mail  Stop  P942 

P.O.  Box  1663 

Los  Alamos,  NM  87545 

Los  Alamos  National  Laboratorv 

ATTN:  Dr.  S.  J.  Czuchlcwski 

Mail  Stop  E543 

P.O.  Box  1663 

Los  Alamos.  NM  87545 

Los  Alamos  National  Laboratorv 

AT  TN:  Dr.  J.  Mack 

Mail  Stop  L554 

P.O.  Box  1663 

Los  Alamos.  NM  87545 


Mission  Research  Corporation 
A  1  I  N:  Dr.  1).  Mitrovhch 
Plasma  Sciences  Division 
I  _2o  Randolph  Road.  SL 
.  Vbbuqucrque.  NM  S"lO(» 
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19.  Naval  Research  Laboratory 
ATTN:  Dr.  J.R.  Greig  (4763) 

4555  Overlook  Avenue.  SW 
Washington.  DC  20375 

20.  Naval  Research  Laboratorv 
ATTN:  Dr.  R.  Hubbard  (4790) 
4555  Overlook  Avenue.  SW 
Washington,  DC  20375 

21.  Naval  Research  Laboratorv 
ATTN:  Dr.  A.  Wahab  Ali  (4700.1) 
4555  Overlook  Avenue,  SW 
Washington,  DC  20375 

22.  Naval  Research  Laboratory 
ATTN:  Dr.  R.  Pechacek  (4760) 
4555  Overlook  Avenue.  SW 
Washington.  DC  20375 

23.  Naval  Research  Laboratory 
ATTN:  Dr.  D.  Murphy  (4760) 
4555  Overlook  Avenue.  SW 
Washington,  DC  20375 


4.  Naval  Research  Laboratorv 
ATTN:  Dr.  M.  Raleigh  (4760) 
4555  Overlook  Avenue.  SW 
Washington,  DC  20375 

5.  Naval  Research  Laboratorv 
ATTN:  Dr.  R.  Fernsler  (4770) 
4555  Overlook  Avenue.  SW 
Washington,  DC  20375 

6.  Naval  Research  Laboratorv 
ATTN:  Dr.  B.  Ilui  (4790)  ' 
4555  Overlook  Avenue,  SW 
Washington.  DC  20375 

7.  Naval  Research  Laboratory 
ATTN:  Dr.  G.  Joyce  (4790) 
4555  Overlook  Avenue.  SW 
Washington.  DC  20375 

2S.  Naval  Surface  Weapons  Center 
White  Oak  Laboratorv 


Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
ATTN:  Dr.  A.  Smith  (H23) 
10901  New  Hampshire  Avenue 
Silver  Spring.  MD  20903-5000 

Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
ATTN:  Dr.  II. C.  Chen  (R41) 
10901  New  Hampshire  Avenue 
Silver  Spring,  MD  20903-5000 

Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
ATTN:  Dr.  H.S.  Uhm(R41) 
10901  New  Hampshire  Avenue 
Silver  Spring.  MD  20903-5000 

Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
ATTN:  Dr.  R.  Fiorito  (R41) 
10901  New  Hampshire  Avenue 
Silver  Springs,  MD  20903-5000 

Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
ATTN:  Dr.  J.  Smith  (R41) 

1<)901  New  Hampshire  Avenue 
Silver  Springs.  MD  209"3->ooo 

Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
ATI  N:  Dr.  D.  Rule  (R41) 
l09o]  New  Hampshire  Avenue 
Silver  Spring.  MD  20903-5000 

Naval  Suriace  Weapons  Center 
White  Oak  Laboratory 
ATTN:  Dr.  M.J.  Rhce  (R41) 
l09oj  New  Hampshire  Avenue 
Silver  Spring.  ML)  2<)9o3-5<»oo 

Oflice  ol"  Naval  Research 
A  I  I  N:  CDR  James  Olfutt 
l'»30  Last  Green  Street 
Pasadena.  C.\  9 1  loo 

O! iice  of  Naval  Research 
A  I  I  N:  (  DR  R.  Swallord 
soi)  \.  ()uinc\  Street 
Arhncton.  V.\  2221" 


mm 


38.  Admiral  R.L.  Topping 

Space  and  Naval  Warfare  Svstcms  Command 
SPA  WAR-06 

Washington.  DC  20363-5100 

39.  IX DR  E.  Turner 
p  MW- 145 

SPA  WAR 

Washington,  DC  20363-5100 

40.  Sandia  National  Laboratories 
ATTN:  Dr.  C.  Ekdahl  (1272) 

P.O.  Box  5800 
Albequerque,  NM  87 185 

41.  Sandia  National  Laboratories 
ATTN:  Dr.  R.  Lipinski  (1272) 

P.O.  Box  5800 
Albequerque.  NM  871S5 

42.  Sandia  National  Laboratories 
ATTN:  Dr.  M.  Mazarakis  { 1272) 

P.O.  Box  5800 
Albequerque,  NM  871 S5 

43.  Sandia  National  Laboratories 
ATTN:  Dr.  D.  Hasti 

P.O.  Box  5800 
Albequerque,  NM  87185 

44.  Sandia  National  Laboratories 
ATTN:  Dr.  C.  Lrost 

P.O.  Box  5800 
Albequerque,  NM  S71S5 

45.  Sandia  National  Laboratories 
ATTN:  Dr.  J.  Freeman  (1241) 

P.O.  Box  5800 
Albequerque,  NM  87185 

46.  Sandia  National  Laboratories 
ATTN:  Dr.  G.T.  Leifeste  (1241) 

P.O.  Box  5801) 

Albequerque.  NM  S 7 1 S5 

47.  Dr.  KAY.  Struve 

Lawrence  Livermore  National  Laboratorv 
P.O.  Box  Sos 
Livermore.  GA  lM55* » 
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